Evaporative Cooling of a Photon Fluid to Quantum Degeneracy 
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We demonstrate that the process of evaporative cooling, as associated with the cooling of atomic 
gases, can also be employed to condense a system of photons giving rise to coherent properties of 
the light. The system we study consists of photons in a high-quality Fabry- Perot cavity with photon 
interactions mediated by a nonlinear atomic medium. We predict a macroscopic occupation of the 
lowest energy mode and evaluate the conditions for realizing a narrow spectral width indicative of 
a long coherence time for the field. 



The phenomenon of coherence has played a crucial role 
in many areas of physics. An extraordinarily long coher- 
ence time is the fundamental property that distinguishes 
laser light from ordinary light [l|. The extensive coher- 
ence properties of matter have recently been investigated 
in trapped ultracold atoms by studying de Broglie matter 
waves in Bose-Einstein condensates [2|, y, |4| . Coherence 
is a ubiquitous phenomenon: many other systems, such 
as a collection of polaritons [5| , also show similar coherent 
attributes. 

In the case of atoms, coherence presents itself when 
the atomic gas has been condensed to quantum degen- 
eracy. Condensation is often achieved using the tech- 
nique of evaporative cooling in which the highest energy 
atoms are forced to escape the trap, and those that re- 
main rethermalize, thereby reducing the temperature. In 
a nonequilibrium system, continuous evaporative cooling 
may generate a coherent atom laser [6]. In this arti- 
cle, we demonstrate how evaporative cooling and Bose- 
stimulated emission can be used to condense a photon 
fluid into a quantum degenerate superfluid. 

In a general sense, a photon fluid is a collection of spa- 
tially localized photons which, through their interactions 
via a nonlinear medium, exhibit hydrodynamic or similar 
fluidic behavior. A condensed photon fluid should pos- 
sess superfluid properties such as coherence, phase rigid- 
ity, quantized vortices, a critical Landau velocity, and 
the Bogoliubov dispersion [7|. This is distinct from a 
normal laser, where the effective photon interactions are 
so weak that they are unimportant. Another significant 
difference with usual lasers is that a population inversion 
of an internal atomic state is not necessary to generate 
the coherent light. Instead the phase-space compression 
to produce a macroscopically occupied mode relies on an 
inversion of photon cavity mode populations with respect 
to the Planck distribution. Since the relevant modes of 
the system are determined by the structural properties 
of the cavity and not atomic quantities, the frequency of 
the coherent field could be highly tunable. 

Two high-reflectivity mirrors placed close together 
form a Fabry-Perot cavity and support many photon 
modes. The closeness of the mirrors creates a cavity 
mode volume which has a high aspect ratio, with large 



energy gaps between modes which have adjacent longitu- 
dinal quantum numbers and small energy gaps between 
modes which have adjacent transverse quantum num- 
bers. We utilize this to pump the cavity such that only 
one longitudinal quantum number is relevant, and the 
modes considered are distinguished only by their trans- 
verse degrees of freedom. Photons in vacuum interact 
weakly, therefore, it is necessary to incorporate a nonlin- 
ear medium, such as an atomic Rydberg gas or nonlin- 
ear crystal, into the cavity to allow for atom-mediated 
photon-photon interactions. Since the longitudinal de- 
gree of freedom is frozen out due to a small cavity length, 
an intriguing link exists between this system and conden- 
sation in two dimensions with the well-known physics of 
the Berezinskii-Kosterlitz-Thouless transition, which has 
recently been investigated in atomic gases [8|, |9[ . 

The photon fluid in a Fabry-Perot cavity is governed 
by the Hamiltonian, 
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where hi is the annihilation operator for a photon in mode 
i with energy Tia^, and Tf- kl is the scattering or collision 
(C) rate for photons from modes k and / into modes i and 
j. There are two kinds of interactions, as illustrated in 
Fig. [TJ Terms which do not change populations give rise 
to a self-interaction energy and leave the photons in the 
same mode they entered in, e.g. ajajaiao. Population 
transfer is due to recombination terms, such as o^a^aiCLi- 
Photons can also enter and leave the cavity through 
irreversible pumping and decay, which is not described 
by the Hamiltonian in Eq. (pQ). Pumping and decay can 
be incorporated into the system dynamics through the 
quantum master-equation formalism, 
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where p is the density matrix operator of the system, Tf 
is the pumping rate (P) and Tf is, the decay rate (D) 
for mode z, and the Lindbladt superoperator is 



L[d] = 2dpd ] - 6 ] dp - pd ] 6 . 
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FIG. 1: (Color online) The five primary processes that oc- 
cur in the Fabry- Perot cavity are presented. Photons enter 
the system from a pumping source. While in the cavity, the 
photons recombine and self- inter act. Finally, the photons can 
leave the cavity through evaporation or output coupling. 



The relative pumping rates as a function of the frequency 
of the system modes are dependent on the pumping 
source, which, for simplicity, we take here to be constant 
across all relevant modes, Yf — > Y p '. 

In contrast, evaporative cooling requires a decay rate 
from the cavity which is strongly dependent on the en- 
ergy. A mirror with a frequency dependent reflectivity, 
which allows modes to decay at different rates, could be 
utilized. A reflectivity which drops for large frequen- 
cies can create conditions where high energy photons 
leave the cavity faster than low energy photons. Since 
a high energy photon carries away more than the aver- 
age photon energy, subsequent rethermalization through 
photon interactions should reduce the system tempera- 
ture, in much the same way as in the evaporative cooling 
of atoms. 

We implement energy dependent reflectivity by divid- 
ing the decay rates into two classes. High energy photons 
(H) are described by a decay rate of Y H and low energy 
photons (L) by a decay rate of T L , with Y H > Y L . The 
large decay rate of the high energy modes means they 
can be eliminated in the following way. Consider an in- 
teraction in which photons from modes i and j create a 
photon of low energy in mode k and a photon of high 
energy in mode I. The rapid decay of mode I causes this 
process to be effectively irreversible. This is described 
most simply by extracting such processes which involve 
a high energy mode from the reversible H term in Eq. O 
and adding them back to the density matrix evolution as 
an irreversible term of superoperator form, 
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where Y E = (Y c ) 2 /Y H is the evaporation rate (£"), taken 
to be constant. The system Hilbert space is now reduced 
and the high energy modes no longer appear in the theory 
explicitly. 

Even with these simplifications, the problem is still in- 
tractable in general. The dimensionality of the Hilbert 



space grows exponentially with the number of photons in 
the lower energy modes. Most of this complexity is unin- 
teresting from the point of view of evaporative cooling of 
photons since it arises from all possible pathways to redis- 
tribute and entangle photons amongst the lowest energy 
modes. Consequently, we divide the system into what we 
will refer to as plaquettes, each plaquette consisting of a 
pair of modes i and j, which may interact to produce a 
photon in the lowest energy mode k = and a photon in 
a high energy mode I which is lost through the mirrors 
rapidly. Adding the contributions from D distinct pla- 
quettes as an incoherent sum of pumping rates into the 
ground state, the density matrix equation of motion in 
the interaction picture is, 
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where the sum over (i,j) implies modes i and j are from 
the same plaquette. 

Amongst the populations, the evolution time scales are 
not all equivalent. If only one or no atoms are in a plaque- 
tte, a slow evolution takes place on time scales governed 
by Y p and Y L . If two atoms are in a plaquette, however, 
and providing Y E ^> T p , T L , they will rapidly collide and 
form a ground state photon, and a high energy photon 
which will evaporate away. Consistent with this is that 
we may neglect the possibility of three or more photons 
in the plaquette. 

The basis states needed then describe the number of 
photons in each plaquette and the number of photons in 
the lowest energy mode. The precise occupation of the 
modes is not important, only how many plaquettes have 
a given number of photons. For instance, an arbitrary 
basis state |\£) can be uniquely identified with a reduced 
computational state described by, 
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where n is the number of ground state photons, a is the 
number of plaquettes with no photons, b is the number of 
plaquettes with one photon, and c and d are the number 
of plaquettes with two photons in different modes or the 
same mode, respectively. 

The photon number distribution is an important quan- 
tity to characterize the photon state since a coherent 
state has a Poissonian distribution while a thermal state 
has an exponentially decaying distribution. The popula- 
tions of all states |\£) that can be expressed as |n, abed) 
are given by 



Pn,abcd = (*|p|*)- 



(7) 



We can now take advantage of the natural separation 
of time scales by adiabatically eliminating the fast pro- 
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FIG. 2: (Color online) Photon number probabilities of the 
lowest photon mode of steady states with constant number of 
low energy photons. As the pumping power per mode is de- 
creased and more modes are pumped, denoted by decreasing 
temperature, the distributions change from an exponential 
decay associated with a thermal system to nearly Poissonian 
as associated with a superfluid system. 



cesses. This is implemented by solving for the fast vari- 
ables in steady state and substituting back into the equa- 
tions for the slow variables. The populations of states 
with a plaquette with two photons are then approxi- 
mately given by their steady state values, 
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as determined by Eq. [5] and Eq. The evolution of the 
remaining populations is then determined by, 



dP n ,a 

dt 
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where P n?a = P n , a (D-a)oo- 

The number distribution of photons in the low energy 
mode is highly dependent on the incoherent pumping rate 
and the number of plaquettes pumped. The pumping 
rate can be designated by an effective temperature, T, 
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where n is the number of photons that would be in each 
mode if the modes were isolated and non-interacting and 
uj is the bare mode energy. This temperature corresponds 
to the occupancy of the external pumping modes. In or- 
der to make an effective comparison between different 



temperatures, we hold constant the average number of 
photons in the lowest energy mode, and vary both the 
temperature and, correspondingly, the energy width of 
the pump in order to maintain this condition. It would 
be expected that for an intense, narrow-bandwidth pump 
connected only with the lowest energy mode of the Fabry- 
Perot cavity, the high temperature and incoherent char- 
acteristics of the pump would be imprinted on the pho- 
tons in that mode. However, for a weaker but broader 
pump, in which a wider spectrum of modes are pumped 
but at a slower rate so that the average photon number 
remains the same, the lowest energy mode is not as in- 
fluenced by the temperature of the pumping modes. The 
occupation of the lowest energy mode is due mainly to re- 
combination and stimulated emission which should create 
a coherent system in direct analogy with the evaporative 
cooling of atoms. 

Figure [2] presents the number distributions of photons 
in the lowest energy mode during steady state operation 
for systems with 20 photons in the low energy mode. 
With a large effective temperature, as denoted in the 
figure by /c#T = 21^, the number distribution is es- 
sentially what would be expected for a thermal source. 
As the effective temperature decreases, by decreasing the 
pumping strength and increasing the number of pumped 
plaquettes, the effects of bosonic amplification dramat- 
ically increase. This leads to a nearly Poissonian num- 
ber distribution associated with zero temperature con- 
densates. The key to converting the incoherent thermal 
pump into a coherent system is that the occupation of 
the lowest energy mode comes from stimulated emission 
from higher energy modes. 

The spectrum of the photons is also an important 
quantity to examine as it provides information about the 
coherence time and spectral width. Since the mode sep- 
aration is greater than the collisional interaction energy 
and the system is in steady state, the fluctuation spec- 
trum is determined by, 
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which is the Fourier transform of a two-time correlation 
function. The correlation function can be expressed as, 



<4(*)a (0))= ^v^C n){i) (t), 

n,{i} 



(13) 



(11) with, 



C nAi} {t) = Tr[p\n, {i})(n - 1, {i}\(t)a (0)] , (14) 

where n is the number of ground state photons and 
{%} represents the configuration of the remaining modes. 
From the quantum regression theorem, the equation of 
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FIG. 3: (Color online) The spectrum is presented with van- 
ishing collision rates. As the photon number increases, the 
spectral width narrows. In the inset, the inverse of the spec- 
tral width is proportional to average photon number as ex- 
pected from the Schawlow-Townes linewidth. 



FIG. 4: (Color online) The spectrum is presented for systems 
with different collision rates. Larger collision rates lead to 
wider spectra. The spectrum center is also increased from 
the bare mode frequency due to mean- field interaction shifts. 



motion for the quantity C n ^ is given by. 



dt 
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where the matrix elements are determined from the off- 
diagonal density matrix equation of motion, 



9p n {i},{n-l){i} 
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as given by Eq. [51 As this is a reduced set of the full 
density matrix equations of motion, it is merely necessary 
to replace these matrix elements into Eq. ([T5]) with the 
initial condition C nj ^}(0) = \/nPn,abcd(0), where \n, {i}} 
corresponds to the reduced basis state |n, abed). 

There are two driving forces of decoherence, incoherent 
transfer through the mirrors and phase dispersion due 
to collisions, which both limit the coherence time. The 
decoherence due to transfer through the mirrors can be 
examined by setting a vanishing collision rate. Due to 
bosonic amplification, the photons are able to support 
a state with a spectral width much narrower than that 
given by the linewidth of the photon decay. 

Figure [3] presents the spectrum for vanishing collision 
rates for systems with 20 pumped plaquettes and dif- 
ferent average photon numbers. As the photon number 
is increased due to a higher pumping rate, the spectral 
width decreases due to Bosonic amplification. The spec- 
trum can be much narrower than the width associated 
with the decay rate through the mirrors. As shown in 
the inset of Fig. 02 the spectral width, denned as half 
width at half maximum, is inversely proportional to the 
average photon number as expected by the Schawlow- 
Townes linewidth |l|. 

The effects of non- vanishing collision rates are of course 
an important aspect of the system. Figure [4] presents 



the fluctuation spectrum for a system with an average of 
100 photons as the collision rate is increased. There is 
a mean-field interaction shift of the frequency since the 
energy of the photons is no longer just the bare mode 
energy. The spectral width rapidly grows as the collision 
rate becomes larger. One possible way to remedy this 
problem is to create a system that has a small collision 
rate for the low energy photons but still possesses a large 
recombination rate for the high energy photons. Such 
a system can be manufactured by creating a spatially 
dependent nonlinear medium. The low energy photons 
are more localized to the center of the cavity than the 
high energy photons. Therefore, if the medium is less 
dense in the center of the trap than toward the edges, the 
cavity will allow for fast recombinations of high energy 
photons to promote rapid evaporative cooling but small 
collision rates for low energy photons to promote coherent 
evolution. 

The condensation of a photon fluid using evaporative 
cooling provides a new mechanism for creating a high- 
intensity narrow-spectrum coherent beam of light. Unlike 
a laser, the photon fluid, due to the presence of interac- 
tions, is expected to be superfluid and should, for exam- 
ple, exhibit signature properties such as phase rigidity of 
the order parameter, and the ability to support quantized 
vortices. 
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